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Abstract 

We identify three independent two-fold parameter degeneracies (5,0 13 ), 
sgn(<5m 2 ]) and (023,tt/2 — @23) inherent in the usual three-neutrino analy¬ 
sis of long-baseline neutrino experiments, which can lead to as much as an 
eight-fold degeneracy in the determination of the oscillation parameters. We 
discuss the implications these degeneracies have for detecting CP violation 
and present criteria for breaking them. A superbeam facility with a baseline 
at least as long as the distance between Fermilab and Homestake (1290 km) 
and a narrow band beam with energy tuned so that the measurements are per¬ 
formed at the first oscillation peak can resolve all the ambiguities other than 
the (@23 , 7 t /2 — @23) ambiguity (which can be resolved at a neutrino factory) 
and a residual (5 ,ir — 5) ambiguity. However, whether or not CP violation 
occurs in the neutrino sector can be ascertained independently of the latter 
two ambiguities. The (5 , n — 5) ambiguity can be eliminated by performing 
a second measurement to which only the cos 5 terms contribute. The hier¬ 
archy of mass eigenstates can be determined at other oscillation peaks only 
in the most optimistic conditions, making it necessary to use the first oscilla¬ 
tion maximum. We show that the degeneracies may severely compromise the 
ability of the proposed Super JHF-HyperKamiokande experiment to establish 
CP violation. In our calculations we use approximate analytic expressions for 
oscillation probabilitites that agree with numerical solutions with a realistic 
Earth density profile. 


I. INTRODUCTION 


The up/down asymmetry of the neutrino flux (originating from cosmic ray interactions 
with the atmosphere) at SuperKamiokande is now a 10 a effect. A compelling interpretation 
of this result is that neutrinos have mass and oscillate from one flavor to another. The 
atmospheric neutrino deficit is explained as a consequence of z/ /t —> u T oscillations with 
almost maximal amplitude and mass-squared difference, dm 2 l ~ 3x 1CV 3 eV 2 [JIJ. The K 2 K 
experiment || with a baseline of 250 km has preliminary results that are in agreement with 
this interpretation. Oscillations of zy, to v e as an explanation of the atmospheric anomaly 
are ruled out by the CHOOZ |[]| and Palo Verde [[1|] reactor experiments, which place a 
bound on the amplitude smaller than 0.1 at the 95% C.L. in the Smh region of interest. The 
MINOS 0, ICARUS || and OPERA [[7| experiments are expected to come online in 2005 and 
study aspects of the oscillations at the atmospheric scale 0 . The low energy beam at MINOS 
will allow a very accurate determination of the leading oscillation parameters. ICARUS and 
OPERA should provide concrete evidence that ry —> z/ r oscillations are responsible for the 
atmospheric neutrino deficit by identifying tau neutrino events. 

Measurements of electron neutrinos from the Sun also provide strong evidence for neu¬ 
trino oscillations. The flux of electron neutrinos from the Sun observed in several different 
experiments is smaller than the Standard Solar Model 0 (SSM) prediction by a factor of 
1/3-1/2. The recent SNO charged-current measurements show that u, 
explain the v e flux suppression relative to the SSM 
angle (LMA) and small matter effects (Srn^i 


ry iT oscillations 
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The solution with a large mixing 
5 x 10 ” 5 eV 2 and amplitude close to 0.8) has 


emerged as the most likely solution to the solar neutrino problem [11] . This solution will be 
tested decisively by the KamLAND reactor neutrino experiment 
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There are several parameter degeneracies that enter the determination of the neutrino 
mixing matrix which can be removed only with future oscillation studies with superbeams 
or neutrino factories. See Table [I] for a sample of proposed baselines. A notable example 
is the U e 3 (=sin 6 , i 3 e - * 5 ) element. Only an upper bound exists on 0i 3 , nothing is presently 
known about the CP phase 5, and the two always appear in combination in the mixing 
matrix. It is the breaking of such degeneracies that will be of concern to us in this work. 

In Section II we identify all the potential parameter degeneracies in the mixing matrix. 
We restrict our attention to the 3x3 matrix that describes the mixing of active neutrinos, 
setting aside the possibility that the atmospheric, solar and LSND |T3| data may require the 
existence of a fourth neutrino that is sterile. The parameter ambiguities are connected with 


TABLE I. Baseline distances in km for some detector sites (shown in parentheses) for neutrino 


beams from Fermilab, Brookhaven, JHF, and CERN. 


Fermilab 

Beam source 

Brookhaven 

JHF 

CERN 

150 (Frejus) 


350 (Cornell) 

295 (Super-K) 


730 (Soudan) 



730 (Gran Sasso) 

1290 (Homestake) 


1200 (Seoul) 


1770 (Carlsbad) 

1720 (Soudan) 



2100 (Beijing) 

2640 (San Jacinto) 

2540 (Homestake) 



2900 (SLAC) 

2920 (Carlsbad) 
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not only neutrino mixing but also the neutrino mass pattern; we pay particular attention to 
the implication of these ambiguities for the detection of CP violation. In Section III, within 
the context of a superbeam experiment |14|, we present methods by which all but one of 


these degeneracies can be resolved, and argue that the remaining ambiguity can be settled at 
a neutrino factory [[L5|] . We also discuss the implications of the degeneracies on the proposed 
Super JHF-HyperKamiokande experiment HI , which would have a 4 MW proton driver and 
a 1 Mt water cerenkov detector (40 times larger than SuperKamiokande). We summarize 
our results in Section IV. In an appendix we provide a complete set of approximate analytic 
expressions for the oscillation probabilities that are useful for superbeams and neutrino 
factories, and define their domain of validity by making comparisons with numerical solutions 
of the evolution equations. 


II. PARAMETER DEGENERACIES 

In this section we identify the three types of parameter degeneracies that can occur 
in the three-neutrino framework when —> v e and > v e oscillation probabilities are 
used to extract the neutrino parameters. We use approximate formulas [|T7],|T^] for neutrino 
propagation in matter of constant density to illustrate the degeneracies. In each case we 
discuss the implications for detecting CP violation. 


A. Oscillation probabilities in matter 


The neutrino flavor eigenstates u a (a = e, /i, r) are related to the mass eigenstates 
Vj (j = 1,2,3) in vacuum by 


v* = J2 Kj u i > 

j 


( 1 ) 


where U is a 
described by 


unitary 3x3 mixing matrix. The propagation of neutrinos through 
the evolution equation flKiyiZU 




.dv a 





rn„ 


+ 


A 


2 E,, 2 E, 


^ae^0e I ^0 


( 2 ) 


where x = ct and A/2E V is the amplitude for coherent forward charged-current v e scattering 
on electrons, 


A = 2V2 G f N e E v = 1.52 x 10~ 4 eV 2 V e p ( g/cm 3 )E u (GeV) , (3) 

Here N e is the electron number density, which is the product of the electron fraction Y e (x) 
and matter density p(x). In the Earth’s crust and mantle the average matter density is 
typically 3-5 g/cm 3 and Y e ~ 0.5. The propagation equations can be re-expressed in terms 
of mass-squared differences 
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where 8m 2 k = m 2 — rri 2 .. The neutrino mixing matrix U can be specihed by 3 mixing angles 
(023,012)013) and a CP-violating phase 8. We adopt the parameterization 


U 


* C13C12 C13S12 Si3e lS \ 

-C23S12 - S 13 S 23 C 12 P S C23C12 - Si 3 S23Sl2e l<S C13S23 

V S23S12 - Si 3 C 23 Ci2e* <5 -S23C12 - Si 3 C23Si2e i<5 C13C23 J 


(5) 


where Cjk = cos 6jk and Sjk = sin 6 ^. In the most general U, the % are restricted to the 
first quadrant, 0 < < 7 t/ 2 , with 5 in the range 0 < 5 < 2i r. We assume that v 3 is 

the neutrino eigenstate that is separated from the other two, and that the sign of dm 2 


31 


can be either positive or negative, corresponding to the case where v 3 is either above or 
below, respectively, the other two mass eigenstates. The magnitude of 8rn 2 M determines 
the oscillation length of atmospheric neutrinos, while the magnitude of 8rn 2 n determines 
the oscillation length of solar neutrinos, and thus IhmlJ -C If we accept the likely 

> 0 and we can restrict 0 12 to the 


conclusion that the solar solution is LMA [O 


then 8m 21 


range [0, 7 t/ 4 ]. It is known from reactor neutrino data that #13 is small, with sin 2 20 13 < 0.1 
at the 95% C.L. [[|. Thus a set of parameters that unambiguously spans the space is 8m 31 
(magnitude and sign), 8m | l5 sin 2 20 i 2 , sin 6 * 23 , and sin 2 26 * 13 ; only the 0 23 angle can be below 
or above 7 t/ 4 . 

I 11 the context of three-neutrino models the usual method proposed for detecting CP 
violation in long-baseline experiments with a conventional neutrino beam is to measure the 
oscillation channels —> v e and —> P e (or v e —>• and P e —>■ for a neutrino factory). 
Both leading and subleading oscillation contributions must be involved and the oscillations 


must be non-averaging for CP-violation effects | 2 T|] . For illustrative purposes we use the 


constant density matter approximation, although in an exact study variations of the density 
along the neutrino path should be implemented. Approximate formulas for the oscillation 
probabilities in matter of constant density in the limit -C A, already exist in 


the literature [ p~7| , p~8| . We adopt the form in Ref. |T 8 [], where 6*13 is also treated as a small 
parameter and the mixing angles in matter are found in terms of an expansion in the small 
parameters 6* 13 and 8m\ x . We introduce the notation 


A = \8m 2 31 \L/AE u = 1.27\8m 2 31 /eV 2 \{L/km)/ (P^/GeV), 
A = | A/8m 2 1 


311 


a = | Sm^/Sm 


311 


( 6 ) 

(7) 

( 8 ) 


Up to second order in a and $ 13 , the oscillation probabilities for 8m\ x > 0 and 8'w 2 2] > 0 are 


P(u M —> u e ) = x 2 f 2 + 2xyfg(cos8 cos A — sin 8 sin A) + y 2 g 2 , (9) 

P(P M —>■ u e ) = x 2 f 2 + 2xyfg(cos 8 cos A + sin 8 sin A) + y 2 g 2 , (10) 

respectively, where 

x = sin 6*23 sin 26 13 , ( 11 ) 

y = a cos 6*23 sin 29 12 , (12) 

/, / = sin((l A) A )/(1 =p A), (13) 

g = sin(AA)/A. (14) 
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The coefficients / and / differ due to matter effects (A^O). To obtain the probabilites for 
Srri 2 j < 0, the transformations A —> —A, y —> —y and A —> —A (implying / <->■ —/ and 
(7 —> —g) can be applied to the probabilities in Eqs. H) and (|T 0 |) to give 

P(zy t —> i/ e ) = x 2 ] 2 — 2xyfg{cos5cos A + sin A sin A) + y 2 g 2 , (15) 

P(h /t —> P e ) = x 2 / 2 — 2xyfg{cos S cos A — sin <5 sin A) + y 2 g 2 . (16) 

For a T-reversed channel, the corresponding probabilities are found by changing the sign 
of the sin 5 term. In Eqs. ( 0 ), ( 0 ), ( 0 ), and (JTTJ) we have assumed 8 rn 2 n > 0 , which is 
what one expects for the LMA solar solution; for < 0, the corresponding formulae are 
obtained by y —> —y. These expressions are accurate as long as $i 3 is not too large, and 

they are valid at E u > 0.5 GeV (A > 0.04(3 x 10 ” 3 eV 2 /|Am 2 1 |)) for all values of dm^ 

currently favored by solar neutrino experiments. We expand on the domain of validity of 
these equations in the Appendix. The corresponding expansion in a and 6 13 in a vacuum 
can be found by the substitutions f,f,g—> sin A. 

For reference, the conversion from A and A to L and E v is shown in Fig. [Tj. For neutrinos 
with dm 2 ± > 0 or anti-neutrinos with dm^ < 0, A — 1 corresponds to an MSW resonance. 
For neutrinos, it can be shown that the choice A = 1/2 maximizes both the sin A and cos A 
terms for a given A; for anti-neutrinos the A that maximizes the sin A and cos 5 terms varies 
with A. 

We make two observations regarding the approximate probability formulas above, the 
consequences of which are discussed below: 

(i) Both terms that depend on the CP phase 5 vanish when <7 = 0 , i.e., at AA = nn, 
where n is an integer. The y 2 term also vanishes in this case, so that only the x 2 term 
survives. 

(ii) The cos A term vanishes when A = (n — -) 7 r, while the sin 5 term vanishes when 
A = 7171. 


The above statements are true for both neutrinos and anti-neutrinos. 

The first observation implies that there is no sensitivity to the CP-violating phase 5 
when ( N e )L = / N e dL is an integer multiple of \[2ii/Gp, where (N e ) is the average value of 
N e for the neutrino path. Numerically, for n — 1, this condition is 

(N e )L ~ 16275 km , (17) 

or, for the Earth’s density profile, 

P~ 7600 km. (18) 


This distance has a simple physical interpretation: it is the characteristic oscillation wave¬ 
length due to the matter interaction [|19|]. Furthermore, the condition in Eq. ([l7j) is inde¬ 
pendent of all oscillation parameters. It is also independent of E u . It has often been noted 
that CP violation is strongly suppressed in long baseline experiments of order 7300 km 
(nominally the distance from Fermilab to Gran Sasso); we see that this is a universal effect 
that occurs because L is close to the oscillation length due to matter. Furthermore, the term 
proportional to y 2 also vanishes, which means that there is also no dependence on Srri^ or 612 
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0 71/4 71/2 371/4 71 


A 

FIG. 1. Contours of (a) L and (b) E u versus A and A, for = 3 x 10 -3 eV 2 . 

at this distance, at least to second order in the small parameters. Therefore this distance is 
especially well-suited for measuring $13 without the complications of disentangling it from <5, 
$ 12 , or dmh. For baselines greater than about 4000 km the constant density approximation 
loses accuracy (see results in the Appendix), so that the critical distance in Eq. ( |T8D is not 
exact, but does explain semi-quantitatively the weakness of CP violating effects near that 
distance. 

The second observation relates to the relative strength of the sin<5 and cosh terms in 
P(/y /t —> u e ). In short L , low E u experiments the matter effects are small and the leading 
terms of the oscillation probability are given by the vacuum formulas. Then L and E v 
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can be chosen such that only the explicitly C'P-violating sin 8 term survives (e.g., when 
A = 7 r/ 2 ), and CP violation can be measured directly by comparing P(z/ M —> u e ) and 
P(P M —> P e ) (although even for L ~ few 100 km there are small matter corrections that 
must be considered). However, as is evident from Eqs. (||) and HD, when # 13 is small the 
relative strengths of the sin 5 and cos 6 terms in the presence of large matter corrections at 
longer L can be selected by an appropriate choice of A in exactly the same way as in the 
short L , vacuum-like case. That is, the 5 dependence with matter effects included can be 
made pure sin <5 for 

/q y in-3 eV 2 \ 

L/E u ~ (2 n - 1)(410 km/GeV ) -—- , (19) 

V r m 3i\ ) 

where n is an integer Q The only caveat is that matter corrections are much larger for 
longer L and the accuracy of the determination of 5 may be more subject to knowledge of 
the electron density. 


B. Orbits in probability space 

We assume that sin 2 2 # 23 and are well-determined (perhaps at the few percent 

level or better) by a survival or —» v T measurement ||, and that #12 and hm 2 ] are also 
well-determined (KamLAND should be able to measure the parameters of the solar LMA 
solution to the few percent level [[1^]). Then the remaining parameters to be determined are 
8, #i 3 , and the sign of Sm^ (the sign of Sm^ is positive for LMA). 

The usual proposal for testing CP violation in the neutrino sector is to measure both 
zy t —* v e and —* u e probabilities. As 6 varies for given #i 3 and sgn(5m 2 1 ), an elliptical 
orbit will be traced in P—P space | [H)j . fT?|| . The shape of the ellipse is determined by the 
relative phases of the terms involving 5. We identify three possible cases: 

(i) A 7 ^ mi/2. In this case, both the sin 6 and cos 6 terms are nonzero and the orbit for 
fixed #13 is an ellipse. Each value of 5 gives a distinct point in P—P space for a given 
#i 3 . For A = (n — |)|vr the ellipse has the maximum “fatness” [^], i.e., it is as close 
as possible to a circle given the values of / and /. 

(ii) A = (n — |)7r, where n is an integer. In this case the cos S term vanishes and the 
orbit ellipse collapses to a line. If / ~ / (such as at short L where matter effects are 
small), CP violation is measured directly by comparing the v and D event rates (after 
correcting for the differences in the cross sections and initial flux normalization). 

(hi) A = wr. In this case the sin 6 term vanishes, the ellipse collapses to a line, and CP 
violation is measured indirectly by parametrically determining the value of 8 and not 
by the measurement of a C'P-violating quantity. 


*The misconception that the cos 8 term dominates at large L and E v comes from extending the 
large E u approximation beyond its range of validity, as discussed in Ref. ||l 8 | . 
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There will be two ellipses for each 9 13 , one for each sign of Sm 3l ; they both fall into the 
same class, i.e., if the ellipse for Smh > 0 is Case (ii), the ellipse for Sm 3l < 0 will also be 
Case (ii). 


C. CP degeneracy: ( 6 , @13) ambiguity 


In many cases the parameters (h, 6 * 13 ) can give the same probabilities as another pair of 
parameters (S', 9' 13 ), for fixed values of the other oscillation parameters; this is known as the 
“(<S> # 13 ) ambiguity” |23|. Using Eqs. (|9|) and (|10|), the general formulas for the parameters 
(x', S') that give the same P and P as ( x, 5) for A 7 ^ mr/2 (Case (ii)) are 


/ cV 

X COSO 

x ‘' sin S' 


x cos S + 
x sin 5 — 


if + f)( X 2 ~ x ' 2 ) 

Ayg cos A 
(f - f)(x 2 - x' 2 ) 
Ayg sin A 


( 20 ) 

( 21 ) 


Equations 


x 


/2 


and (| 2 l| ) can be used to derive 

Ayg sin 2A yg sin 2A + xf sin( A — S) + xf sin(A + S) 


x 2 = 


/ 2 + / 2 - 2 //cos 2 A 


( 22 ) 


from which S' can then be determined from Eq. (;20i) or m- In particular, a set of parameters 
which violates CP (sin S' 7 ^ 0) can be degenerate with another set of parameters, with a 
different 6 * 13 , that conserves CP (sin S = 0). It can be shown that in all cases real solutions 
exist for x', so there will be an ambiguity between two sets of oscillation parameters if 
| sin A' | < 1. Therefore we conclude that the use of a monoenergetic beam at a fixed L 
will necessarily entail parameter ambiguities if only the channels —> v e and —> v e are 
measured. An example is shown in Fig. ^a for A = ‘Air/A. 

When A 7 ^ mr/2, the (S,9i 3 ) ambiguity can give a degeneracy between CP violating 
(CPV) and CP conserving (CPC) solutions. If sin <5 = 0 in Eq. (^l[), then sinh' is not zero 
if / P /; hlie difference can be large if / and / differ substantially due to large matter effects. 
For example, in Fig. |2]a the prediction for (P, P) for (sin 2 2013, <5) = (0.01, 0) is identical to 
that for (0.00298, 1.48 tt). 

For the Cases (ii) and (iii) above, the ellipse collapses to a line. The ambiguities then 
reduce to x' = x (see Eq. (|22|) ) and sinh' = sinh in Case (ii) and cosh' = cosh in Case (iii). 
Thus the ambiguity no longer involves 6*13 (and hence in principle d 13 is determined, at least 
as far as the (h, 6 >i 3 ) ambiguity is concerned), but instead is a (S,tt — h) ambiguity (which 
does not mix CPC and CPV solutions) in Case (ii) and a (h, 2ir — h) ambiguity in Case (iii). 
Examples for Cases (ii) and (iii) are shown in Figs. §b and ^]c, respectively. For A ~ mr/2 
the orbit ellipse is very skinny and the ambiguous 9 13 values are close to each other, which 
qualitatively is similar to either Case (ii) or Case (iii). 

Note that in both Figs. 0 b and i the orbit line has a negative slope. It can be shown 
that for A < 1 (i.e., density less than the critical density for resonance) the orbit lines in 
(P, P) space have negative slope for A = mr/2. For A > 1 , the orbit lines for A = mr/2 
have positive slope. 


8 








L = 1290 km, Sm^ = 3x10 -3 eV 2 , Sm^ = 5x10“ 5 eV 2 



FIG. 2. Orbit ellipses showing (<5, ^ 13 ) ambiguity for L = 1290 km with (a) E u = 2.09 GeV 
(A = 37 t/4 ), (b) E u = 3.13 GeV (A = 7 r/ 2 ), and (c) E u = 1.565 GeV (A = n), for sin 2 2 #i 3 = 0.01 
and 0.00298. The other parameters are Srn^ = 3 x 10 -3 eV 2 , Srri^i = 5 x 10 -5 eV 2 , sin 2 2023 = 1 5 
and sin 2 2012 = 0.8. The value of 5 varies around the ellipse. In (b) and (c) the ellipse collapses to 
a line and the ambiguity reduces to a (d, 7 r — 6 ) or ( 5 , 27t — 5) ambiguity, respectively, and different 
values of $13 do not overlap (for the same sgn(5m 2 3 )). 
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D. Mass hierarchy degeneracy: sgn(5m 31 ) ambiguity 


In addition to the (5, 0 13 ) ambiguity discussed above for a given sgn(5m 31 ), in some cases 
there are also parameters (5', d[ 3 ) with 5m 31 < 0 that give the same P and P as (5, #i 3 ) with 
5m 31 > 0. Three examples of the sgn(5m 31 ) ambiguity are shown in Fig Furthermore, 
there is also a (5' ambiguity for Sm 31 < 0, so in principle there can be a four-fold 
ambiguity, i.e., four sets of 5 and 9 13 (two for 5m 2 1 > 0 and two for Srri 2 { < 0) that give the 
same P and P. 

As with the (5, 0 13 ) ambiguity, the sgn(5m| 1 ) ambiguity can mix CP conserving and CP 
violating solutions. For example, in Fig. |3|a the prediction for (P,P) for (sin 2 29 13 , 5) = 
(0.01,0) with 5m 31 = 3 x 10~ 3 eV 2 is identical to that for (0.0138, An/3) with Srri 2 ] = 
—3 x 10 -3 eV 2 . 

Although the general equations for the sgn(5m 31 ) ambiguity are somewhat messy, for 
the case A = (n — \)i r the values of {xf S') for Sm 31 < 0 that give the same P and P as 
(x, 5) for Sm 31 > 0 are determined by 


,j 2 _ x 2 (f 2 + P - ff) — 2yg(f — /)x sin 5 sin A 

OC — 


x! sin 5’ = x sin 5 


ff 

P + P-ff x 2 p + pf-f 


ff 


sin A ff 2 yg 


If sin 5 = 0 then Eq. 


reduces to 


sin S' = — x 


ff 


P + P f-f _ 

ff 2yg S mA\p + p-ff' 


(23) 

(24) 


(25) 


which is not zero if / 7 ^ /, i.e., whenever there are matter effects, so there is a potential 
CPC/CPV confusion as long as the right-hand side of Eq. (^) has magnitude less than 
unity. It is possible to have S' = 7t/ 2 when 5 = 0, i.e., CPC can be confused with maximal 
CPV (see Fig. |b). 

The ambiguity between parameters with 5m 31 > 0 and 5m 31 < 0 occurs only for some 
values of 5, and does not occur at all if matter effects are large enough (i.e., L and 9 13 are 
large enough) [|4|,^]. The conditions for the existence of this ambiguity will be discussed 


further in Sec. |III Bj . Note, however, that the sgn(5m| 1 ) ambiguity can still confuse different 
values of 5 and 6* 13 even for A = mr/2 (see, e.g., Figs. |3|b and|3]c), unlike the (5, 9\ 3 ) ambiguity 
where 9 33 is removed from the ambiguity for A = n 7 r/ 2 . 


E. Atmospheric angle degeneracy: (6*23 , vr/2 — 023 ) ambiguity 

There is yet another ambiguity in the determination of 5 and 9 1 3 , which involves the 
value of 623. In practice it is only sin 2 29 23 that is determined by a survival measurement 
(for now we ignore matter corrections to which are relatively small for oscillations 

involving active flavors), so 623 cannot be distinguished from 7 t /2 — 6*23. The effect of this 
degeneracy can be seen by interchanging sin 6*23 and cos 6*23 in Eqs. (|TT| ) and (JT^). For 
623 — vr/4 (the favored solution from atmospheric data) the ambiguity vanishes, but for 


10 














L = 730 km, Sm^ = 5x10 5 eV 2 



FIG. 3. Sgn^TOy]) ambiguity for L = 730 km with (a) E u = 3.54 GeV (A = 7r/4), (b) 
E v = 1.77 GeV (A = vr/2), and (c) E v = 0.885 GeV (A = n). The other parameters are 
^ m 2 i = 5 x 10 -5 eV 2 , sin 2 2023 = 1, and sin 2 2012 = 0.8. 
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sin 2 2023 — 0.9 it can have a sizable effect, since in this case sin 2 023 = 0.35 and cos 2 023 = 
0.65. Three examples of the 0 2 3 ambiguity are shown in Fig. |4|. The 0 2 3 ambiguity can 
also mix CPC and CPV solutions; for example, in Fig. Eb the prediction for (P, P) for 
(sin 2 20 13 , sin 0 2 3 , 5) = (0.01,0.585,0) is identical to that for (0.00107, 0.811, 47 t/ 3 ). 

As with the sgn(5m 2 1 ) ambiguity, the equations for the 023 ambiguity are rather messy 
in the general case. For the special case A = (n — |) 7 r, we have 


sin 2 20 (o = sin 2 2013 tan 2 0 2 3 + 


7 13 


a 2 g 2 sin 2 20 12 

77 


(1 - tan 2 0 23 ), 


sin 20(3 sin S' = sin 20 13 sin S + 


ag(f ~ f) sin 2012 cot 20 23 


// 


sin A 


(26) 

(27) 


where ( 5 , 0i 3 ) are the parameters that give a certain (P, P) for 0 < 0 23 < 7 t/ 4 and (S', 9[ 3 ) are 
the parameters that give the same (P, P) for 7r/2 — 0 23 . We see that even for A = nir/2 the 
023 ambiguity can mix CPC and CPV solutions, since sin S = 0 does not necessarily imply 
sin S' = 0. Furthermore, even for A = nn/2 the 623 ambiguity mixes solutions with different 
0i3 (see Eqs. (|^) and (^7|), and Figs. [i]b and £|c), unlike the (5,013) ambiguity where 0 i3 is 
removed from the ambiguity for A = rm/2. 

Since a; is a small parameter (and possibly even small compared to sin 20 i3 ), the numerical 
uncertainty in S due to the 623 ambiguity is generally small, of order 0.077T or less, when 
A — (n — |)7r, L = 2900 km, Srn^ = 10 -4 eV 2 , and sin 2 20 i 3 = 0.01. This effect is of order 
the expected experimental uncertainty in S lpTi| . |2?i|| . The size of the S ambiguity decreases 
with decreasing matter effect (smaller L) and with decreasing Sm^i, so for a wide range 
of parameters the CPC /CPV confusion from the 623 ambiguity is not too severe. On the 
other hand, the sin 2 20 13 confusion is approximately a factor tan 2 0 23 (see Eq. (|25|)), which 
lies roughly in the range | to 2 for sin 2 20 23 > 0.9. 


F. Comments on parameter degeneracies 

In the preceding three sections we have shown that in principle there can be as much as an 
eight-fold ambiguity in determining 5 and 613 from P(iy t —* u e ) and P(z? M —> D e ) at a single L 
and E U1 which comes from the presence of three independent two-fold ambiguities: (5, 0 i3 ), 
sgn( 5 m 3 1 ), and (0 23 ,tt/2 — 0 23 ). For each type of ambiguity there is the possibility of being 
unable to distinguish between CP violating and CP conserving parameters. Measurements 
at multiple L and E u can be used to help discriminate the different degenerate solutions, but 
that would involve extra detectors or a much longer total running time, and probably reduced 
statistics for each ( L , E v ) combination. In the next section we will explore what L and E v 
values do best at resolving these potential degeneracies without resorting to measurements 
at different L and/or E u . We then will examine what L and E u for a second measurement 
can remove the remaining degeneracies. 

In the Appendix we demonstrate that the analytic expressions are accurate for E u > 0.5 
GeV for baselines up to 4000-5000 km. For much lower E v (as low as 0.05 GeV) they are still 
accurate at shorter distances (L < 350 km) if a and 0 i3 are not too large (see the discussion 
in the Appendix). Therefore we expect that the qualitative aspects of the three ambiguities 
are unchanged for short L , low E u experiments such as CERN-Frejus. 


12 









L = 1290 km, Sm^ = 3x10“ 3 eV 2 , Sm^ = 5x10“ 5 eV : 



FIG. 4. (023,tt/2 — # 23 ) ambiguity for L = 1290 km with (a) E u 
(b) E v = 3.13 GeV (A = n/2), and (c) E u = 1.565 GeV (A = ir). 
$ m 2 i = 5 x 10 -5 eV 2 , and sin 2 29\2 = 0.8. 


= 2.09 GeV (A = 3vr/4), 
The other parameters are 
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TABLE II. Possible neutrino beam energies E v (in GeV) versus baseline (in km) and A that 
will convert the (5, $ 13 ) ambiguity to a simple (< 5 , 7 r — 5) ambiguity, for 5m\i = 3 x 10 -3 eV 2 . For 
other values of 5mf n , E v scales proportionately with <Only values of E u > 0.5 GeV are 
considered. 


A 

300 km 

730 km 

1290 km 

1770 km 

2100 km 

2600 km 

2900 km 

7r 

2 

0.73 

1.77 

3.13 

4.29 

5.12 

6.34 

7.03 

3n 

2 


0.59 

1.04 

1.43 

1.71 

2.11 

2.34 

57r 

2 



0.63 

0.86 

1.02 

1.27 

1.41 

77T 

2 




0.61 

0.73 

0.91 

1.00 

9n 

2 





0.57 

0.70 

0.78 


III. RESOLVING PARAMETER DEGENERACIES 


A. Resolving the {5,6 13 ) ambiguity 


As discussed in Secs. [11 B| and m the choice A = nn/2 causes the orbit ellipse in 
(P, P) space to collapse to a line and the {5, # 13 ) ambiguity reduces to one involving only 5, 
i.e., the combination of P and P gives a unique value of $13 (at least for one sign of Am 2 ,). 
Furthermore, since <5 only becomes confused with ir — 5 (in Case (ii) of Sec. [11 Q or 27r —<5 (in 
Case (hi)), CP conserving solutions never become mixed with CP violating ones. Case (ii) 
(with A — (n— tj) 7 t) has another advantage in that the —> v r oscillation is approximately 
maximal (see Eq. (|A3|) ), which would facilitate a better measurement of sin 2 29-23 and hm 2 
Therefore the choice A = (n — \)-k is the best for resolving the (5, 613 


23 a,nu. uno 31 . 

ambiguity. Some 


representative beam energies for particular baselines are given in Table ||. 


B. Resolving the sgn(<5m 31 ) ambiguity 


The parameter degeneracy associated with the sign of 5m 31 can be overcome if there is 
a large matter effect that splits P and P, e.g., if L is sufficiently long and # 13 is not too 
small [^|,^5|]. To determine the minimum value of 9 ± 3 that avoids the sgn(<5m 2 1 ) ambiguity, 
we must first End the region in (P, P) space covered by each sgn(<5m| 1 ), and then determine 
the condition on $13 that ensures the regions for different sgn(<5m| 1 ) do not overlap. 

The orbit ellipse for a given sgn(<5m| 1 ) moves as #13 changes, sweeping out a region in 
(P, P) space. All points on each orbit ellipse (for a given # 13 ) that lie inside the region 
will overlap an orbit ellipse for a different d 13 (this is what leads to the (<5, d 13 ) ambiguity). 
However, the points on the orbit ellipse that lie on the boundaries of the region do not have 
a {5,9 13 ) ambiguity, i.e., there are unique values of 9 13 and 5 for that point. This implies 
that for points on the boundary of the region, x = x' and 5 = 5' in Eqs. (|20|) (p2|) . For 
A 7 ^ 7 m/ 2 , the condition becomes xf sin(A — 5)+xf sin(A + 5) + yg sin 2A = 0. Solving for 
5 and substituting into Eqs. (|9|) and ([L0|) we find the coordinates of the > 0 envelope 
in (P, P) space are given by 


P 


x 


7 


2 + y 2 g 2 


2 y 2 g 2 p sin 2 2A ± 2 ygf^fz{f cos 2A - /) 
f 2 + f 2 — 2f f cos 2 A 


(28) 
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where 


2 = x 2 (f + J 2 - 2//cos2A) - y 2 g 2 sin 2 2A , 


(29) 


and P is found by interchanging / / and letting g —> —g. For (5 m 2 , < 0, the values of 

P and P on the envelope can be found by interchanging P and P. Although the general 
solution is complicated, for the special case A = (n — |) 7 r it is not hard to show that the 
two sgn(<5m| 1 ) regions do not overlap if 


x > 


2yg 

f-r 


(30) 


Note that matter effects split / and /, which decreases the minimum value of value of x 
(and hence of sin 2 26 * 13 ) needed to avoid any sgn^m^) ambiguity. Because x oc sin 26* 13 and 
y oc ( 5 / 77 . 2 !, the corresponding minimum value of sin 2 26* 13 increases as the square of Srn^i- 
The minimum values of sin 2 2# 13 /((5m | 1 ) 2 are plotted versus E v for various values of L in 
Fig. |. For A = 7 t/ 2, the sgn(<5m§!) ambiguity would be resolved for sin 2 26* 13 > 0.01 (0.04) 
at L — 1290 km, the distance from Fcrmilab to Homestake, if Smh = 5 x 10 ” 5 (10~ 4 ) eV 2 . 
For L ~ 2600 km (Brookhaven-Homestake or Fermilab-San Jacinto), sgn(<5m| 1 ) can be 
determined for values of sin 2 26* 13 as low as 0.002 (0.008) for 8 m\ x = 5 x 10 -5 (10~ 4 ) eV 2 . 



FIG. 5. Minimum value of sin 2 2613 /that avoids the sgn(<5m 2 !) ambiguity, plotted ver¬ 
sus L for A = 7 t/ 2 (solid curve) and 2>tt/2 (dashed), with = 3 x 10~ 3 eV 2 . The corresponding 
values of E u are marked on the curves. 


Figure 0 shows that A = 3tt/2 would be unsatisfactory in distinguishing sgn(<5m 31 ); 
in fact, measurements at A = (n — |) 7 r do increasingly worse as n increases, as can be 
shown using Eq. ([131). For A = (n — we have |///| = (1 + A)/{l — A)] since A 
is proportional to E v , and E u decreases with n for fixed L , larger values of n will have 
smaller A. Thus, the values of / and / will be closer for larger n, reducing the size of 
the matter effect (at least as far as splitting P and P is concerned). For A = 3tt/2 and 
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8m 2 2 = 5 x 10 -5 eV 2 , the value of sin 2 20i 3 must be greater than about 0.25 which is excluded 
by CHOOZ |J. Even the most optimistic case for A = 37T/2 (which occurs for the highest 
value of <5/77,2!(— 10~ 4 eV 2 ) allowed in the LMA region) requires sin 2 20 13 > 0.06. Thus, 
the proposal of Ref. [|27| to perform experiments at higher n suffers from an inability to 
determine sgn(<5m 31 ). Practically speaking, only n — 1 will provide sufficient discrimination 
for sgn(<5m| 1 ) if A is restricted to the values (n — |)7r. We henceforth restrict ourselves to 
this case. 

By combining A = 7 t/ 2 with a sufficiently long L, the combined four-fold ambiguity 
involving <5, ffi 3 , and sgn^mgj can be reduced to a simple (<5, 7r — <5) ambiguity that in 
principle determines whether CP is conserved or violated. Some possibilities are shown 
in Figs. | and 0. As Fig. || shows, for 8m\ x = 5 x 1CT 5 eV 2 , L = 1290 km is (barely) 
sufficient to distinguish sgn(h7?7| 1 ) for sin 2 20i 3 as low as 0.01. However, for 8m 21 = 10 -4 eV 2 , 
L > 2000 km is needed. In practice, experimental uncertainties and uncertainties on the 

increase the likelihood of having a sgn(<5r77 2 1 ) ambiguity, so that 


matter distribution 28 


a separation of the two regions greater than the size of the experimental uncertainties is 
required. 


C. Resolving the (023,7 t /2 — $23) ambiguity 


Even if A is chosen to mitigate the effects of the (<5, ffi 3 ) ambiguity, and L is chosen 
long enough to eliminate the sgn(<5m| 1 ) ambiguity, there still remains the (P 23 , 7t/2 — 6 * 23 ) 
ambiguity. If 023 — 7t/4 this ambiguity disappears, and choosing L and E u such that A = ir/2 
leaves a simple (<5, n — 8) ambiguity. If 0 2 3 deviates from 7 t/ 4, then there does not appear to 
be a judicious choice of a single L and E u that can resolve the 0 23 ambiguity. 

The problem in resolving the 023 ambiguity lies in the fact that in the leading term in 
P(/y /i —> v e ) and —> P e ), sin20i 3 is always paired with sin0 23 (see Eqs. @ and (|T0|)), 
and so if there are two values of 023 derived from the measured value of sin 2 202 3 , there will 
be two corresponding values of sin 2 20 13 . Since sin 2 2 0 23 can be as low as 0.9, the two values 
of sin 2 0 23 can be as far apart as 0.35 and 0.65, and therefore the ambiguity in sin 2 20 13 can 
be as large as a factor 1.86 at leading order (Eq. (p6| ) in the limit that a is small). The 
next-to-leading term in the probabilities in Eqs. (||) and (|10]) is proportional to sin20 23 , 
and therefore cannot resolve the ambiguity. The last term in P(u M —> u e ) is proportional 
to cos 2 02 3 , so that the relative weighting of the last term compared to the leading term is 
affected by the value of sin0 23 . However, the last term is suppressed by a 2 , and is generally 
much smaller than the leading term (at least for sin 2 20 13 > 0.01, the approximate region 
where superbeam experiments will be able to probe). Hence, even measurements at a second 
L and E v would likely be unable to resolve the 0 23 ambiguity if it exists (i.e., 623 not close 
to 7 t/4). 

If one could also measure P(i/ e —> u T ) (see Appendix A for an approximate analytic 
expression), then a comparison with P(u e —> u^) should determine whether 0 23 is above 
or below 7t/4; the leading term in P(zz e —> v T ) can be obtained from the leading term in 
P(v, —> v e ) by the replacement of sin0 23 by cos02 3 . A v e —> v T measurement could be done 
in a neutrino factory; in fact, a neutrino factory may be the only practical way to resolve 
the 023 ambiguity, if it exists. A neutrino factory experiment also provides energy spectrum 
information that could be helpful in resolving parameter ambiguities . 
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(a) L= 1290 km 


E v = 3.13 GeV 
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5m 31 - 5 = 71/2 
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FIG. 6 . Resolution of combined (<5, $ 13 ) and sgn^m^) ambiguities when A = 7 r/ 2 , for (a) 
L = 1290 km, (b) L = 1770 km, and (c) L = 2900 km, with \8rnf iX = 3 x 10 -3 eV 2 , 
|<5m 2 i| = 5 x 10 -5 eV 2 , and sin 2 2023 = 1- 
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FIG. 7. Same as Fig. |(| except for = 10 4 eV 2 . 
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D. Measurements at a second L and/or E v 


As we have demonstrated, measurements at a single L and E u cannot resolve all param¬ 
eter ambiguities. A second experiment at a different L and/or E v , with a different value of 
A, is required for this purpose. The best sets of L and E v are those that are complementary, 
i.e., the second experiment should provide the clearest distinction between the parameter 
ambiguities of the first experiment. In this section we discuss three possible scenarios, each 
with measurements at two L and E v combinations. 


1. Scenario A 

In this scenario, the first measurement would be done at Ai = 7 t/ 2 (with L/E v given 
by Eq. (|1^)). As discussed earlier, this choice isolates the sin 5 term, removes 6 *i 3 from the 
# 13 ) ambiguity, and the remaining (5, 7 r — 6) ambiguity does not mix CPC and CPV 
solutions. These L/E u values also give a large disappearance, which facilitates the pre¬ 
cision measurement of <5m 3 j and sin 2 26 * 23 . The baseline L should be large enough to avoid 
the sgn(5m| 1 ) ambiguity (L > 2000 km assures this for sin 2 26*i 3 > 0.01). Representative 
values of L and E u are given in Table |TJ. Measuring P(y,j —> v e ) and P(P^ —> u e ) at one 
such L and E u should determine sgn^m^), 6 *i 3 (modulo the 6*23 ambiguity, if present), and 
whether or not CP is violated (as discussed in Sec. |II E| . the existence of a 6* 23 ambiguity 
will not give a large amount of CPC /CPV confusion). 

The second measurement should be one that best resolves the (S,7r — 6) ambiguity. In 
principle, A 2 = 7 r, which eliminates the sin 8 terms in the probabilities and leaves only cos 8 
terms, gives the maximal separation of <5 and 7 t — 5. Thus, the first measurement gives sin <5, 
the second gives cos 8, from which the value of 5 may be inferred. Furthermore, if 6 * 1 3 is 
determined from the first measurement, then both P and P would not have to be measured 
in the second measurement; one is sufficient to determine 5. Whether one used neutrinos 
or antineutrinos in the second measurement would be determined by which gave the larger 
event rate, taking into account neutrino fluxes, cross sections, and oscillation probabiltities. 
If 8rn\ [ > 0, then neutrinos would be best for the second measurement, due to the larger 
flux and cross section; for <5m 31 < 0 , antineutrinos may be the better choice if the matter 
enhancement is enough to overcome the lower flux and cross section for antineutrinos. If 
both the first and second measurements are done at the same L, then A 2 = 7r means that 
the appropriate energy in the second experiment is E% = E\j 2. 

In practice, there are other values A 2 that are not close to 7 t/2 that could potentially 
work for the second measurement. The optimal A 2 also depends on the particular values of 
/, / and g at the various L and E v , as well as on neutrino parameters that are currently 
unknown (Sm^, 6 * 12 , and $i 3 ). We do not pursue the optimization here. 


2. Scenario B 

If the first measurement is done at an L that is not large enough to resolve the sgn^m^) 
ambiguity, then the second measurement must be tailored to both determine sgn(dm 2 j) and 
resolve the (5,tt — S ) ambiguity, i.e., it must break a four-fold degeneracy. As discussed 
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TABLE III. Possible sets of neutrino beam energies and baselines that will resolve the four-fold 
parameter ambiguity when the measurements are done at shorter L (such that the sgn^mjjj) 


ambiguity is not resolved in either experiment). 


Fixed L (km) 

Ei (GeV) 

E 2 (GeV) 

300 

0.73 

0.355, 0.375 

730 

1.77 

0.835, 0.940 

Fixed E (GeV) 

Li (km) 

L 2 (km) 

0.73 

300 

575, 630 

1.77 

730 

1295, 1700 


above, the choice A 2 = vr determines 5, but it can be shown that at shorter L an approximate 
degeneracy with parameters of the opposite sgn(<5m| 1 ) remains (e.g., see Fig. [I|c). Another 
example is shown in Fig. |^a, where the near degeneracy of parameters with the opposite 
sgn(5m| 1 ) remains for some values of <5 (while the crosses in Fig. |8|a are well-separated in 
the second measurement, the boxes are not). However, at A 2 = 7r/(l ± A), either / or / 
vanishes (depending on the sign of 5m 31 ), and the four ambiguous solutions occupy four 
separate regions in (P, P) space, as shown in Fig. |8|b. Although the four regions in Fig. |]b 
overlap somewhat, when the point of one degenerate solution is in the overlap region the 
points of the other three degenerate solutions are not (the crosses and boxes are always well- 
separated). Thus the four-fold ambiguity involving (5,7r — 5) and sgn(<5m 31 ) will always be 
resolved. A disadvantage of Scenario B is that because either / or / is zero in the second 
measurement, P and P tend to be smaller, so that event rates may be somewhat lower than 
for other values of A. 

Some examples are given in Table |T|. For instance, if the first baseline is L\ = 730 km 
(Fermilab to Soudan), then A x = 7 t/ 2 for E = 1.77 GeV, and two possibilities for a second 
experiment with the same beam energy are L 2 = 1295 and 1700 km, which serendipitously 
are very close to the distances from Fermilab to Homestake and from Fermilab to Carlsbad 
(or Brookhaven to Soudan). 

In practice, narrow band beams are not monoenergetic. However, values of A 2 close to 
7T/ (1 ± A) also give reasonably good separation of the ambiguities, as long as A 2 is not close 
to 7r. If the fractional beam spread is more than \A\, a slightly different average value of A 2 
might be preferable, to ensure that no significant part of the beam has A 2 too close to 7 r. 

Fig. | summarizes the possibilities for Scenerios A and B, showing P„ versus L for the 
first measurement done at Ai = 7t/2 (solid curve) and a possible second measurement at 
A 2 = 7T (Scenario A, dotted curve) or 7r/(l ± A) (Scenario B, dashed curves). 


3. Scenario C 


This scenario uses the fact that the probabilities are insensitive to the parameters of 
the 5m 2l scale at L ~ 7600 km, as noted in Sec. |IFBj. If the first measurement of P and 


P were done at L ~ 7600 km, d 13 would be determined (modulo the 8 2:i ambiguity), and 
because the distance is large enough sgn^m^) would also be determined from the large 
matter effect. A second measurement of P and P could then be done at an L and E v such 
that A = (2 n — 1)7t/ 4, which gives the maximum “fatness” of the orbit ellipse j[2^ 


see 


Sec. |IIB|), which in turn should best distinguish different values of 5. One disadvantage of 
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L = 730 km , E, =1.77GeV 




FIG. 8. Values of P and P in a second measurement when there is a four-fold degener¬ 
acy in the first measurement for (a) E 2 = 1.18 GeV (A 2 = 37r/4) and (b) E 2 = 0.94 GeV 
(A 2 = 7r/(l + A) = 0.947t), with L\ = L 2 = 730 km and E\ = 1.77 GeV (Ai = 7r/2). Each 
curve (solid, dashed, dotted, dash-dotted) represents one of the four solutions that are degenerate. 
Points labelled by the same symbol (crosses or boxes) correspond to solutions that are degenerate 
with each other in the measurement at L\ and E\. 

Scenarios B and C compared to Scenario A is that both P and P must be determined in 
both measurements. 


4■ Discussion of scenarios 

Although the three scenarios discussed above are not necessarily the only solutions to 
the ambiguities, in each case one measurement is chosen to eliminate one or more of the 
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FIG. 9. Values of L and E u for a first measurment at Ai = 7 r /2 (solid curve) and a second 
measurement at A 2 = vr (Scenario A, dotted) or 7 r/(l ± A ) (Scenario B, dashed), which breaks the 
parameter degeneracy in each case. 

parameters from the ambiguities, leaving the second measurement to resolve only the re¬ 
maining ambiguities. In this sense, they are cleaner measurements. Scenario A would appear 
to be more favorable, since in principle the first measurement alone could determine sin <5, 
sgn^mjj), and 0 13 (modulo the 6*23 ambiguity), and thus determine whether or not CP is 
violated. Also, as discussed in Sec. [II E| . even if there is a 6*23 ambiguity, the magnitude 
of the CPC/CPV confusion appears to be relatively small for the usual range of neutrino 
parameters considered. 


5. Implication for JHF experiments 


The proposed SuperJHF-HyperKamiokande experiment [TTj| satisfies the requirements 
for the first experiment of Scenario B. The plan is to have a neutrino energy such that 
A is at the first peak of the oscillation for L = 300 km; if A is not exactly on the peak 
(e.g., if Srrif^ = 3 x 10 ~ 3 eV 2 ), a long narrow ellipse results instead of a straight line (see 
Fig. PH|) . Because the distance is relatively short, the sgn(Srrif n ) ambiguity is not likely to 
be resolved since there is considerable overlap of the two sgn(5m 31 ) ellipses (see Fig. |IO|b). 
For example, for sin 2 2ffi 3 = 0.01 the point for 5 = 0 with Sm^ > 0 is nearly the same as 
the point for 6 = 1 . 187 T for 8 rri\ x < 0 (Eq. |25], which measures the size of the CPC/CPV 
confusion for the sgn(hm 2 ,) ambiguity, gives the about the same numerical result). The 
expected 90% C.L. uncertainty in 8 in this case is about 0.077T near 5 = 0 [|16j, so we 
see that the sgn(<5m| 1 ) ambiguity caused by the matter effect would seriously impede a 
proper measurement of 6 , although there is the possibility that the SuperJHF-HyperK 
experiment might measure a point (P, P) that was sufficiently outside the overlap region, 
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thereby determining sgn(5m| 1 ) [| 22 |. A possible ( 6 * 23 , 7 r /2 — 623 ) ambiguity also remains, 


which could lead to a corresponding ambiguity in 0 13 , as shown in Fig. Wf- 

Even though JHF may not sit exactly on the peak of the oscillation (i.e., A = 7 r/ 2 ), 
Fig. [l(] shows that the three ambiguities discussed here are present. Also, Fig. [T 0 | a shows 
that the ambiguity in 6 13 is relatively small (of order 10% or less) if A is close, but not 
exactly equal, to tt/2. Thus as long as L/E v is chosen so that the oscillation is close to the 
first peak, we expect that the scenarios discussed here for determining the neutrino mass 
and mixing parameters will be valid. 


IV. SUMMARY 

There is an eight-fold, (5, 6 * 13 )-sgn(< 5 r? 7 | 1 )-( 6 l 2 3 ,vr/2 — Q 23 ) degeneracy affecting the neu¬ 
trino mixing matrix determined in long-baseline neutrino oscillation experiments. If sin 2 20 23 
is almost unity as is favored by current Super-K and K2K data, this is reduced to a four-fold 
ambiguity. To break this four-fold ambiguity to a simple (£, tt — 5) ambiguity which does not 
interfere with a determination of whether or not CP is violated in the neutrino sector, we 
find that an experiment should be performed at the first oscillation maximum corresponding 
to A = 7 t/ 2 and at a baseline of at least about 1300-2000 km, depending on the value of 
Sm^. Representative values of L and E v that yield A = 7 t/ 2 are given in Table |Id|. 

The obvious advantages of choosing A = 7 t /2 are that z/ M —> v e transitions are nearly 
maximal even when matter effects are accounted for and the v M —> u T oscillation (which has 
small matter effects) is maximal, allowing a precise measurement of sin 2 20 23 and 5m 31 . By 
choosing A = 7 t/ 2 , the (5, 6 * 13 ) degeneracy represented by the P-P ellipse collapses to a line 
leaving a (8,7: — 5) ambiguity (see Fig. |2j) which unambiguously determines whether or not 
CP is violated. The central reason for the choice of the first oscillation peak over other 
peaks is that practically speaking, the sgn(<5m 2 1 ) ambiguity can be resolved only for this 
peak (see Fig. H). The other peaks succeed in eliminating this ambiguity only for the smallest 
values of 5m 21 in the LMA region and for sin 2 2ffi 3 close to the CHOOZ bound. As shown 
in Figs. H|-[7] to remove this ambiguity simultaneously with the (5, 6 * 13 ) ambiguity requires 
that the baseline be at least 1300 km for sin 2 2ffi 3 > 0.01 and <5m 21 = 5 x 10 -5 eV 2 . For 
lower values of sin 2 26* 13 and/or higher values of longer baselines than this are needed. 

The exciting aspect of an experiment at A = n/2 and a sufficiently long baseline is that all 
degeneracies other than the ( 6* 23 , 77/2 — 623 ) degeneracy can be broken to a harmless (5, tt — 8) 
ambiguity with only a single baseline and energy. The remaining (8,tt — 8) ambiguity can 
be removed by making a second measurement at A = tt which leaves only cos 8 terms in the 
probabilities and provides the maximal separation between 8 and n — 8. The (0 23 , 7 t /2 — 6 * 23 ) 
degeneracy cannot be eliminated even with measurements at a second baseline and energy 
because in the leading term in P(z/ M —> u e ) and P(P M —>► P e ), sin 26*13 is paired with sin 6* 23 
(see Eqs. (0) and (|T0|)) . Fortunately, the mixing of the CPC and CPV solutions arising 
from this degeneracy are of order or smaller than the experimental uncertainty in 8, thereby 
making it less severe. Only a neutrino factory, which offers the unique ability to compare 
P(v e —> Vp) and P(u e —> u T ), can disentangle sin 2 d 13 from sin 6* 23 and End whether 6* 23 is 
less than or greater than 7 t/4. 

If it is not possible to have an experiment with L sufficiently large to hnd sgn(<5m 31 ), 
a second experiment is necessary to simultaneously resolve the (5, 7 t — 8) ambiguity and 
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L = 300 km , E v = 0.70 GeV 



FIG. 10. Examples of the three types of ambiguities for the proposed SuperJHF-HyperK ex¬ 
periment |L6| with L = 300 km and E u = 0.7 GeV: (a) (5, # 13 ) ambiguity, (b) sgn^m^) ambi¬ 
guity, and (c) (# 23 , tt/2 — # 23 ) ambiguity. In each case Sm^i = 5 x 10~ 5 eV 2 , sin 2 2#23 = 1, and 
sin 2 2#i2 = 0.8, unless otherwise stated in the figure. The circle in (b) indicates the size of the 
expected experimental uncertainties JTT 
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determine sgn(<5m| 1 ). One possibility is to chose A = 7 r, but as suggested by Fig. |3]c, 
a sgn(hm| 1 ) ambiguity may still remain when experimental errors are included. If A = 
7 r/( 1 ± A), either / or / vanishes, depending on sgn(5m| 1 ), and the four-fold degeneracy 
breaks into four separate regions as in Fig. [||b. See Table 0 for some examples of how this 
scenario can be implemented. The proposed SuperJHF-HyperK experiment would satisfy 
the requirements for the first measurement of this type; it has the limitation of a possible 
sgn(<5m| 1 ) confusion that leads to an ambiguity in the value of J, which may compromise its 
ability to unambiguously establish CP violation. 

In Fig. |9]we summarize the baselines and energies for two measurements, one at Ai = ir /2 
and another at either A 2 = 7r (if the first measurement can determine sgn(5m| 1 ) and only 
the (h, 7 r — 5 ) ambiguity needs resolution) or A 2 = 7r/(l ± A) (if the first measurement can 
not be performed at a long enough baseline and the four-fold degeneracy (5,7r —h)-sgn(^m| 1 ) 
needs to be broken). Another possibility is to have one measurement at L ~ 7600 km with 
AA X = 7 r and a second measurement with A 2 ~ (2n — 1)7 t /4. If K2K, MINOS, ICARUS and 
OPERA End that # 2 3 is not very close to 7 t /4, a neutrino factory will be needed to resolve 
the (623 ,7 t /2 -623) ambiguity. 

In our analysis we have assumed that Sm^ and sin 2 2623 are known when the experiments 
described here are done. In fact, they will likely be determined only to 10% or so. Once 
these uncertainties are included the minimum value of L required to resolve the sgn^m^) 
ambiguity, e.g., in Scenario A, could be slightly longer than indicated here. Also, because 
is not precisely known, the average neutrino energy will not necessarily be exactly at 
the peak defined by A = 7 t/2 . However, as our analysis of the proposed SuperJHF-HyperK 
experiment shows, only minimal uncertainties in 5 and sin 2 2613 are introduced by these 
factors, and the three principal ambiguities discussed in this paper will be qualitatively 
unchanged. 
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APPENDIX A: 


We provide a complete set of analytic expressions for the off-diagonal probabilities that 
are valid in the regime |A| > |a|, which roughly translates to E u > 0.5 GeV. The diagonal 
probabilities can easily be found from them. The off-diagonal probabilities are for a normal 
hierarchy (in addition to Eqs. (|9]) and (|IUD) 


P(zz e —> z/ T ) = cot 2 9 2 3X 2 f 2 — 2xyfg(cos5 cos A + sin S sin A) + tan 2 # 23 7/ 2 g 2 , (Al) 

P(u e —> D t ) = cot 2 023 X 2 J 2 — 2xyfg(cos5cosA — sin S sin A) + tan 2 d 23 7/ 2 g 2 , (A2) 


and 
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P(v n —> u T ) = sin 2 26*23 sin 2 A 

+a sin 26*23 sin 2A f -- sin 6*13 sin 26 *i 2 cos 20 2 3 sin A — A cos 2 6*i 2 sin 26*23^ • (A3) 

VI — A / 


For P(z/ M —■► u T ), replace A by —A in Eq. (]A^). Note that P(v M —> z/ T ) is independent of 5 to 
0(a). To obtain the probabilites for an inverted hierarchy, the transformations A —»■ — A , 
a —> —a and A —> —A must be made (implying / —/ and g —*• —<7 in Eqs. (M-A2), and 
for the T-reversed channels the sign of the sin 6 term must be changed. 

We now compare the results of the analytic expressions with the numerical integration 
of the evolution equations of neutrinos through the Earth. We integrate the equations along 
a neutrino path using a Runge-Kutta method. The step size at each point along the path 
is 0.1% of the shortest oscillation wavelength given by the scales dm 2 { and A. We account 
for the dependence of the density on depth by using the Preliminary Reference Earth Model 
(PREM) 0. To calculate the analytic probability, we use the average value of the electron 
density along the neutrino path. We provide some values in Table |TV] for the reader’s 
use; they are not indicative of the precision with which the electron density is known. We 
include subleading 6*13 effects, which however are not relevant for sin 2 2^13 of 0(0.01) or 
smaller. They are of importance at #13 for which the CHOOZ limit sin 2 20 13 <0.1 (at 95% 
C. L.) is saturated. The parameters chosen to make this comparison are 8m^ x = 3.5 x 10~ 3 
eV 2 , (5'/7i|, = 5 x 10~ 5 eV 2 , 6*23 = 7t/4, 6* 12 = vr/6, sin 2 26* 13 = 0.01 and E u = 5 GeV. Thus, 
the ensuing comparison is not affected by dropping subleading 6*13 effects. A normal mass 
hierarchy (Smf n > 0) is assumed. We will comment on the comparison involving an inverted 
mass hierarchy. 

Figure ITT] shows P(u^ —>► v e ) and P(P M —> u e ) versus distance for 5 = 0,7t/2,7t, and 
77 t/ 4. The agreement between the analytic formulae (solid lines) and the numerical results 
(dashed lines) is excellent for distances up to about 4000 km. Beyond that, the overlap 
between the lines degrades and for L > 5000 km, the analytic equation completely breaks 
down. The analytic expression for P(P M —* v e ) works for much longer distances than that 
for P(i/ M —► u e ). Analogously, for the inverted mass hierarchy, P (z/ p —> v e ) is valid to longer 
distances than P(h M —> u e ). 

Reference |T8|] claims good agreement between the analytic and numerical results for 
L even larger than 10000 km when a constant density is assumed for the Earth’s density 
profile. The use of a realistic density profile as in the PREM model shows that the agreement 
deteriorates at much smaller distances. 

For the sake of completeness we display the corresponding comparisons for P(r / e —> v T ), 
P(h e —> P T ) and P(r / At —»► v T ) in Figs. [L2] and 0- The parameter values chosen are the same 
as for Fig. 0. The range of validity of Eqs. (|Al|) and (]A|) is the same as for Eqs. (|) 
and O. However, Eq. (|A3| ) agrees almost exactly with the numerical result for the entire 
range considered. This is because matter effects are very small in comparison to the leading 
contribution. For the same reason, P(P M —*► P r ) is almost identical to P(^ —>> v T )- 

In Fig. [1^, we show how well the analytic probabilities P(^ M —*► u e ) and P(P^ —*► u e ) 
agree with the numerical integration for L = 2900 km (the longest baseline emphasized in 
this work) as a function of neutrino energy. The oscillation parameters used are the same as 
for Fig. 0. The precision is remarkable for the spectrum of energies of interest. For shorter 
baselines, the agreement gets even better. 
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We now make some cautionary remarks. Our comparisons were made for a = 0.0143, the 
parameter in which the series was expanded, and sin 2 2# 13 = 0.01 which is assumed to be no 
greater than of O(a). These values are motivated by the existing reactor bounds and global 
fits to the atmospheric and solar data. However, as either of these parameters gets larger, 
the agreement between the analytic equations and the numerical results deteriorates at long 
baselines even if subleading #i 3 effects are included. Conversely, the agreement improves 
with smaller values of a and #13. As a rule of thumb, we recommend that the constant 
density approximation to the probabilities be used only for distances less than 4000 km. As 
can be seen from Fig. for L < 4000 km, the density profile is nearly constant for most of 
the neutrino path, thereby satisfying the implicit assumption (of a constant density profile) 
under which analytic probabilities are valid. 

We have stated that the analytic expressions are accurate for E v > 0.5 GeV for baselines 
of 4000-5000 km. This robust bound can be relaxed for L < 350 km to E u as low as 0.05 
GeV. However, for such low values of E v , the sensitivity of the analytic probabilities to a 
and #13 is high and care must be taken in their use. For example, a comparison with a 
numerical integration is desirable if a and #13 are relatively large. 
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Baseline (km) 

(N e ) 

0-800 

1.284400 

900 

1.320264 

1000 

1.356722 

1100 

1.374506 

1200 

1.474359 

1300 

1.516450 

1400 

1.543904 

1500 

1.563676 

1600 

1.578441 

1700 

1.590218 

1800 

1.599484 

1900 

1.607012 

2000 

1.613276 

2100 

1.618752 

2200 

1.622966 

2300 

1.626689 

2400 

1.630101 

2500 

1.633021 

2600 

1.635155 

2700 

1.636799 

2800 

1.638902 

2900 

1.640515 

3000 

1.641341 

3100 

1.642923 

3200 

1.643719 

3300 

1.644499 

3400 

1.653481 

3500 

1.659367 

3600 

1.663974 

3700 

1.668762 

3800 

1.672956 

3900 

1.676851 

4000 

1.681093 


TABLE IV. The average values of the electron density ( N e ) at baselines for which the analytic 
approximations of the probabilities accurately represent the numerical integration of the evolution 
equations. 
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FIG. 11. P(yn — > v e ) and P{pn &e) versus L for 5 = 0,tt/2,tt,77t/4:. The agreement between 
the analytic formulae (solid lines) and the numerical results (dashed lines) is excellent for distances 
up to about 4000 km. The parameters chosen to make this comparison are 5m\i = 3.5 x 10~ 3 eV 2 , 
<5m|i = 5 x 10“ 5 eV 2 , 023 = vr/4, d \2 = vr/6, sin 2 2013 = 0.01 and E u = 5 GeV. 
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FIG. 12. P{v e —> v T ) and P(^ e —> j/ t ) versus L for the same set of parameters as in Fig. 
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L (km) 

FIG. 13. P(z/^ —> v T ) versus L. The analytic expression (solid line) and the numerical calcula¬ 
tion (dashed line) agree almost exactly for the entire range in L. P {> u T ) is almost identical 
to P(vn —> z'r) because of insignificant matter effects. 
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FIG. 14. P(Vfj, —> t'e) and P(z/„ —> P e ) vs Iq, for T = 2900 km. The oscillation parameters are 
the same as in Fig. [llj. The solid lines (analytic equations) and dashed lines (numerical evaluation) 
are almost undistinguishable. 
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FIG. 15. Density profiles along a selection of of chords of length L passing through the Earth; 
the horizontal axis is the fraction of the total path length. 
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